Motional Narrowing of Electron Spin Resonance in Single- Walled Carbon Nanotubes 
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We have measured electron spin resonance in bulk single-walled carbon nanotubes. The linewidth 
was observed to narrow by as much as ~50%, as the temperature was increased from 3 to 300 K. 
From the linewidth change versus temperature, we find the hopping frequency to be 285 GHz. 
Although the spins are mobile, they are still a finite-level-spin species, as indicated by the Curie law 
temperature dependence of the spin susceptibility. We show that the Dysonian lineshape asymmetry 
follows a three-dimensional variable-range hopping behavior from 3 to 20 K and extract a localization 
length of the hopping spins to be ~100 nm. 

PACS numbers: 76.30.-v, 72.20.Ee, 73.63.Fg 
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Understanding spin dynamics is key to a broad range 
of modern problems in condensed matter physics [iHy] 
and applied sciences 0, [I]- Spin transport is a sensi- 
tive probe of many-body correlations as well as an indis- 
pensable process in spintronic devices. Confined spins, 
particularly those in one dimension (ID), are predicted 
to show strong correlations @, @| and long coherence 
times Single- walled carbon nanotubes (SWCNTs) are 
ideal materials for studying ID spin physics due to their 
long mean free paths and weak spin-orbit coupling [lo| . 
Exotic spin properties in metallic SWCNTs at low tem- 
peratures and high magnetic fields have been predicted, 
including the appearance of a peak splitting in the spin 
energy density spectrum, which can be used to probe 
spin-charge separation in Luttinger- liquid theory [J-|6| . 

One method for studying spin dynamics is electron spin 
resonance (ESR), which provides information on spin- 
orbit coupling, phase relaxation time, spin susceptibility, 
and spin diffusion. Many theoretical and experimental 
ESR studies of SWCNTs have been performed over the 
past decade 0,0, El "12 1- Unfortunately, experimental 
ESR studies of SWCNTs have been difficult due to sam- 
ple issues, such as the presence of trace amounts of fer- 
romagnetic catalyst particles and SWCNT sample inho- 
mogeneity. Substantial confiicts have thus emerged in 
the literature, such as the temperature (T) dependence 
of the spin susceptibility and ESR lineshape |llMl7|. 



Here, we show that the probed spins in SWCNTs are 
localized species that undergo phonon-assisted hopping. 
Specifically, wc find that the ESR linewidth becomes 
smaller as T is increased, a hallmark of the phenomenon 
of motional narrowing fisl [l9| . This spin mobility ac- 
counts for the Dysonian lineshape [2^ seen throughout 
the full T range examined. In addition, we observe that 
the spin susceptibility displays a Curie law (oc 1/T) be- 
havior, suggesting that, while spins may move about the 
sample, they are still a finite-level system. The Dysonian 
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FIG. 1: (color online) ESR scans of SWCNT powder sample, 
(a) A raw spectrum at 5.0 K across a 6000 G applied field 
range shows a sharp ESR peak on top of a large FMR back- 
ground, (b) Raw detailed scans with Dysonian lineshape fits 
for all temperatures investigated. Each trace is intentionally 
offset for clarity, (c) Comparison of highest and lowest tem- 
perature lineforms. The linewidth. A, B, and a (see text for 
definitions) are given for both traces for comparison. 



lineshape asymmetry parameter, which is proportional 
to the conductance of the probed spins, is shown to fol- 
low the 3D variable-range hopping (VRH) trend. From 
our data, we find the density of states at the Fermi en- 
ergy to be ~ 10^^ states/eV-cm'^, and we estimate the 
localization length of the spins to be ~ 100 nm. 

Data was taken using a Bruker EMX spectrometer us- 
ing a cavity in the TE102 mode in the X-band (9.6 GHz) 
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regime. For temperatures below 125 K, we used 200 /iW 
of microwave power, while for temperatures 125 K and 
above, we used 1 mW. In all cases, we were in the lin- 
ear absorption regime as shown by microwave power de- 
pendent scans performed on the sample across the full 
temperature range. Our sample consisted of acid puri- 
fied laser oven SWCNTs in powder form [2lj. We com- 
pacted the powder by ultracentrifugation after immers- 
ing the SWCNTs in water. This procedure immobilized 
the nanotubes into a condensed pellet without the need 
for direct mechanical compaction. In addition, adsorbed 
gases and water were removed by annealing the sample 
for 24 hours at 115°C under vacuum with a small argon 
gas flow. To protect the sample from contact with the 
outside environment, the 630 ^g (0.24 g/cm^) SWCNT 
sample (-27% of which is SWCNTs equivalent to 
170 /ig or 0.064 g/cm^) was submerged in mineral oil 
under helium gas in a sealed quartz tube [22 1. 

ESR spectra (Fig. [1]) were taken at discrete tem- 
peratures from 3.4 to 300 K. We used 2,2-diphenyl-l- 
picrylhydrazyl (DPPH) as a frequency standard to cal- 
ibrate the applied magnetic field, Hq, and a 1 mM 
CuS04-5H20 standard as a concentration standard to ex- 
tract the spin susceptibility. As seen in Fig.UJa), a broad 
ferromagnetic resonance (FMR) dominates the spectrum. 
The FMR is tentatively assigned to the nickel and cobalt 
catalyst particles remaining in the acid-purified sample. 
By closely scanning Hq around the SWCNT ESR peak, 
a detailed study of the SWCNT ESR lineshape can be 
performed as a function of T, as seen in Fig. [TJb). Both 
the linewidth and peak-to-peak amplitude become larger 
as T is decreased, while the line-center position (g-factor, 
or simply, g) shows little T dependence, staying close to 
the free-electron g-factor of 2.0023 for each different tem- 
perature examined. The ESR line is asymmetric, as seen 
in Fig. [T](c) , having what is often referred to as a Dyso- 
nian lineshape 20] indicating movement of the electrons 
in and out of the Hi perturbing magnetic field. 

To gain further quantitative understanding, we numer- 
ically fit each ESR scan. The slowly varying FMR back- 
ground was fit by a combination of two large-linewidth 
(- 1000 G) Lorentzians. The SWCNT ESR feature was 
fit using the weak form of the Dysonian lineshape 3, 2 
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where y = , is the resonance field, Ai/ is the 

half- width and is equal to , with 7 = and /ib the 
Bohr magneton. The weak form of Dysonian can be used 
here, because the conductivity and diffusion of the elec- 
trons in the SWCNT powder are both low as compared 
to a traditional metal. Nevertheless, unlike traditional 
magnetic resonance where the signal depends entirely on 
the imaginary part of the ac spin susceptibility, x"i the 
Dysonian lineshape is also influenced by the real compo- 
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FIG. 2: (Color online) ESR linewidth as a function of tem- 
perature (T). The steep rise in AH as T is decreased shows 
that the spin dephasing time becomes shorter as the lattice 
energy is decreased. A theoretical fit based on variable range 
hopping, given by the black line, well-describes the data. 



nent, x' ■ Taking the ac susceptibility, x^ to be 

X = x"cosLp + x'sini^. 



(2) 



we define a = tant^, which is a dimensionless measure 
of the relative contribution of the real part (= for tra- 
ditional ESR). This parameter is also a measure of the 
asymmetry of the lineshape, and A/B ss 1 -f a when 
a ^ 1, where A and B are defined in Fig. [Ijc). Fur- 
thermore, it can be related to the electrical conductivity 
of the probed spins, cTspin, as a cx (Xspin, allowing us to 
compare transport and spin resonance results [231] . 

From the numerical fitting, we extract AH, a, and 
the mass spin susceptibility, Xg. Unlike previous stud- 
ies of ESR in nanotubes, AH changes substantially with 
T. As Fig. [2] shows, as T is increased, AH rapidly de- 
creases until —25 K, whereupon the linewidth begins to 
decrease more slowly up to 300 K. This decrease of the 
ESR linewidth with increasing T is consistent with the 
phenomenon of motional narrowing [isl . [Toj . which oc- 
curs because the dephasing time of the spins can change 
as their translational energy is altered. At high T, the 
spins move rapidly, allowing for less time around dephas- 
ing centers, thus reducing the interaction between the 
probed spins and the dephasing centers. This decreased 
interaction gives a longer spin dephasing time (T2), which 
in turn narrows the lineshape. Conversely, at low T, the 
spins are moving more slowly, which broadens the line. 

To understand the observation of motional narrowing 
more quantitatively, we start with a generalized narrow- 
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ing model [2J, [23 : 



Ai7 



(3) 



where AHp is the amphtude of the perturbations and 
AiJe describes the rate of spin motion. Equation ^ 
was originally derived to describe exchange narrowing or 
motional narrowing from spin diffusion. However, spin 
diffusion can be described in terms of phonon-activated, 
variable range hopping with a probability, dvrh, that is 
proportional to exp (-2i?/^ - /\E/k^T) [261^. Here, 
R is the hopping distance, ^ is the localization length, IS.E 
is the average spacing between energy levels, and Ub is the 
Boltzmann constant. Combining this hopping conduc- 
tion with Eq. dSj and adding an offset, Ai/o, gives t29,'33| 
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where AHq is the high-T ( "metallic" ) asymptotic limit of 
the linewidth and A is independent of T. As Fig.[2]shows, 
Eq. (j4]) fits very well to the observed linewidth. We find 
A to be 11.6±0.8 meV-G and AHq = 5.10±0.08 G, with 
the latter quantity giving a T2 spin dephasing time of 
11.1 ns. The activation energy, AE, is 1.18±0.09 meV 
(13.7 K or 285 GHz). From the value of AE, we can 
get an estimate of how much time (on average) each spin 
spends at each hopping location, r = = 558 fs [3ll |. 
If we phenomcnologically take T2 = nr, where n is the 
number of jumps before phase coherence is lost, then we 
can estimate n to be on the order of 10^ hops, where we 
have taken T2 to be ~10 ns. 

The temperature dependence of our linewidth and the 
fitting to Eq. ^ suggests that we are observing ESR 
of localized carriers, Jn a similar fashion to P-doped 
Si 



32, mi and Ge m 



|. In those systems, mo- 
tional narrowing of ESR is seen at low and intermediate 
donor concentrations. However, if the impurity concen- 
tration is increased, more complicated models would be 
needed [3^ 34 1. Additionally, although motional nar- 



rowing has not been observed in SWCNTs before, it is 
reasonable to expect that type of behavior in a system 
such as nanotubes, since electrons in the basal plane of 
graphite, in addition to ID conducting polymers, also 



show motional narrowing [35l-l3< 



The T dependence of Xg scales with 1/T, indicating 
that the spins are a finite-level system (Fig. [3]). Indeed, 
Xg can be fit very well using a Curie law 
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where fi is the spin magnetic moment and N is the num- 
ber of probed spins per unit volume. The Curie coeffi- 



cient, C 



is found to be l.lliO.Ol x 10"'^ emu- 



K/ g, w hich agrees with the values from other stud- 
ies [l6| . We do not fit the last two points (at 3.7 
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FIG. 3: (Color online) The mass spin susceptibility, Xg, as a 
function of inverse temperature follows a Curie-law down to 
4.0 K, as given by the black best fit line. 



and 3.4 K), since it appears as though Xg becomes con- 
stant below 4.0 K; further work is needed to establish the 
low-T spin susceptibility trend. In contradiction with Pe- 
tit et al. [12|, we find no evidence of Pauli-like behavior; 
if we were observing it, we would expect some sort of 
finite spin susceptibility to remain at high T, as is seen 
in other ID systems like protonated polyaniline 3^ 40| . 
The Curie law behavior of the spin susceptibility has been 
also observed in lightly doped Si and Ge systems 32. 33|. 
as well as intermediately doped Si [4l|. At higher dop- 
ing concentrations, however, deviations from the Curie 
law are seen, with high doping concentrations displaying 
Pauli-like trends [UllHiil. 

Using the spin susceptibility and the measured mass 
of the sample, we can estimate the number of spins per 
nanotube. To simplify calculations, we choose to imagine 
the 170 ^g sample made up of entirely (10,10) SWCNTs. 
Given the Curie constant, we have 3.04 x lO^'^ probed 
spins in the sample or 1.14 x 10^^ spins/cm'^ since our 
sample volume is 2.66 x 10^"^ cm'^. A given 40 carbon 
atom unit ceU for a (10,10) SWCNT is 0.246 nm in length. 
If each nanotube is estimated to have a length of 1 /ini, 
we have 0.58 spins per SWCNT, or a Curie coefficient 
of 3.6 X 10^25 cmu-K/SWCNT. This spin density per 
nanotube suggests that nearly all nanotubes ('^2/3 of the 
total mass) contribute one spin to the observed signal. 

To gain deeper insight into the spin hopping mecha- 
nism, we examined the asymmetry Dysonian lineshape 
parameter, a, which is proportional to the conductance 
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of the probed spins. In particular, we were interested 
to see if a followed a VRH behavior at low T, which is 



mathematically given as [27|, |28 1 



a = ao X exp 



(6) 



where Tq is the characteristic temperature and d is the 
dimensionality of the system. As Fig. 2] shows, ln(a) fol- 
lows a linear trend with T"*, indicating that the spins 
probed by ESR are undergoing phonon-assisted hopping 
in 3D from 3.4 to 20 K. We find from our fit that Tq is 
17.9±5.5 K and that ao is 1.20. The asymptotic limit of 
the a parameter, ao, approaches 1 as T — >• 0. This trend 
is expected, since the asymmetry of the ESR signal is 
caused by thermally-activated hopping: as the phonon 
density decreases, so does the asymmetry of the Dyso- 
nian. The localization length, £, o f the electronic wave- 
function can be found from To 281 
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18.1 



kBToDiE^ 



1/3 



(7) 



where D(Ef) is the density of states around the Fermi en- 
ergy, Ep. We can estimate D{E-p) by treating the defect 
density of states as having an energy separation that can 
be roughly estimated by A£;. Thus, D{Ey) ~ - 
10^^ states/cm'^-eV, where we are utilizing the spin den- 
sity extracted from the Curie constant, N{E-p) = 1.14 
X 10^^ spins/cm'^. From Eq. ([7|), we estimate ^ to be 
'-^lOO nm, similar to previous measurements of defect- 
induced localization lengths in SWCNTs Q ■ The spac- 

ing of defects, R^, can be estimated by (^-/V) , or 
~28 nm. A d — VRH behavior is exp ected in this 
wavefuncation-overlap regime (i?d < £,) [27|. In addition, 
in this regime, exchange effects may also be important, 
but defect concentration dependence is needed to inves- 
tigate this avenue more fully. 

It is important to note that given the difficulty in ex- 
tracting a, the trend that we are observing can be 
considered robust. We also performed traditional con- 
ductance measurements on a similarly prepared sample, 
and although the conductance clearly showed 3D VRH 
behavior, we believe that our ESR and four-point probe 
conductivity measurements are probing different species, 
since the hopping parameters do not agree and the high- 
T trends are different. 

Finally, we note that the origin of the ESR signal is 
from SWCNTs, despite the fact that amorphous carbon, 
graphite, and other inhomogeneities are present in our 
sample. Unlike graphite, our observed ESR signal has a 
much lower ^-factor. In addition, the T-dependence of 
the linewidth does not follow an EUiott-Yafet relaxation 
mechanism, as it does in graphite (35| . Similarly, we can 
discount the signal arising from amorphous carbon; in 
that system, the linewidth variation with T is not con- 
sistent with our observations 451. 
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FIG. 4: (Color online) The natural log of the asymmetry 
Dysonian parameter, a, plotted against the fourth root of 
inverse temperature. The trend of ln(a) follows a 3D VRH 
behavior, as given by the best fit line in black. 



Although it is clear that the ESR signal arises from 
nanotubes, its microscopic origin is not certain. In a 
separate paper j22[, we advocate that n-type defects are 
essential for the SWCNT ESR. We base this conclusion 
on the observation that ESR signal strongly depends on 
the presence of molecular oxygen (a p-type acceptor), 
an effect we find to be compatible with a compensation 
mechanism. This hypothesis is fully consistent with the 
data we present here. ESR-active defect states would 
explain the Curie-law dependence of Xg^ "^^11 as the 
localized, phonon-assisted hopping of the spins and the 
motional narrowing of /S.H . In addition, our results are 
similar to those seen in doped semiconductors. The wide 
variety of SWCNT ESR results over the past two decades 
may be due to the different defect concentrations, which 
would change both the temperature dependence of the 
spin susceptibility and AH(T). Further work needs to 
be done to elucidate and clarify this hypothesis. 

In summary, we have performed temperature- 
dependent ESR measurements on an ensemble of SWC- 
NTs. We observe motional narrowing of the ESR 
linewidth in nanotubes for the first time. The spin sus- 
ceptibility shows a 1/T-type trend, confirming that the 
probed spins are localized species. An estimate of the 
spin density gives 0.58 spins per nanotube. We observe 
that the Dysonian asymmetry parameter, a, follows a 
trend, which suggests 3D variable range hopping 
as the spin transport mechanism. From the extracted 
parameters, we estimate the localization length to be 
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'^lOO nm with a density of states at the Fermi energy 
of - 10^^ states/eV-cm^. 

This work was supported by the DOE/BES through 
Grant No. DEFG02-06ER46308, the NSF through Grant 
Nos. OISE-0530220 and OISE-0968405, the Robert 
A. Welch Foundation through Grant No. C-1509, the 
Air Force Research Laboratories under contract number 
FA8650-05-D-5807, the W. M. Keck Program in Quan- 
tum Materials at Rice University, the Korean Ministry 
of Education, Science and Technology under the World 
Class University Program (R31-2008-10029), and the 
NIH National Heart, Lung, and Blood Institute through 
Grant No. HL095820. We thank Qimiao Si, Adilet Imam- 
bekov, Robert Hauge, and Vladimir Burka for useful dis- 
cussions. 



[1 

[2: 

[4 

[5: 
[6: 
[7: 
[s: 
[9: 

[lo; 

[11 

[12 
[13 
[14 

[15 
[16 



lkono@r ice^ed u( corresponding author. 
Q. Si, Phys. Rev. Lett. 78, 1767 (1997) 



3191 



85, 3464 



ibid. 81, 

(1998). 

L. Balents and R. Egger, Phys. Rev. Lett. 
(2000); Phys. Rev. B 64, 035310 (2001). 
A. A. Kiselev and K. W. Kim, Phys. Rev. B 61, 13115 
(2000). 

S. Rabello and Q. Si, Europhys. Lett. 60, 882 (2002). 

A. De Martino et al, Phys. Rev. Lett. 88, 206402 (2002). 

B. Dora et al, Phys. Rev. Lett. 101, 106408 (2008). 
D. P. DiVincenzo, Science 270, 255 (1995). 

S. A. Wolf et al, Science 294, 1488 (2001). 
T. Giamarchi, Quantum Physics in One Dimension (Ox- 
ford University Press, Oxford, 2004). 
T. Ando, J. Phys. Soc. of Japan 69, 1757 (2000). 
M. Kosaka et al, Chem. Phys. Lett. 233, 47 (1995). 
P. Petit et al, Phys. Rev. B 56, 9275 (1997). 
A. S. Claye et al, Phys. Rev. B 62, R4845 (2000). 
K. Shen, D. L. Tierney and T. Peitrafi, Phys. Rev. B 68, 
165418 (2003). 

J.-P. Salvetat et al, Phys. Rev. B 72, 075440 (2005). 
V. Likodimos et al, Phys. Rev. B 76, 075420 (2007). 



[17 
[18 

[19 

[20 
[21 
[22' 
[23^ 
[24; 

[25' 
[26' 
[27 
[28^ 



[29 
[30| 

[31 

[32' 
[33; 

[34; 

[35 
[36 
[37 

[38' 
[39 

[4o; 

[41 

[42; 

[43' 
[44' 
[45^ 



B. Corzilius et al, Phys. Stat. Sol. (b) 245, 2251 (2008). 
N. Bloembergen, E. M. Purcell, and R. V. Pound, Phys. 
Rev. 73, 679 (1948). 

R. Kubo and K. Tomita, J. Phys. Soc. Japan 9, 888 
(1954). 

F. J. Dyson, Phys. Rev. 98, 349 (1955). 

P. Nikolaev et al, J. Phys. Chem. C 111, 17678 (2007). 

W. D. Rice et al, arXivlll0.5220vl (2011). 

V. Sitaram et al, Phys. Rev. B 72, 035209 (2005). 

P. W. Anderson and P. R. Weiss, Rev. Mod. Phys. 25, 

269 (1953). 

P. W. Anderson, J. Phys. Soc. Japan 9, 316 (1954). 

A. Miller and E. Abrahams, Phys. Rev. 120, 745 (1960). 

N. F. Mott, Phil. Mag. 19, 835 (1969). 

H. Kamimura and H. Aoki, The Physics of Interacting 

Electrons in Disordered Systems (Clarendon Press, New 

York, NY USA, 1989), 1st ed. 

D. K. Wilson, Phys. Rev. 134, A265 (1964). 

K. Morigaki and T. Mitsuma, J. Phys. Soc. Jpn. 20, 62 

(1965). 

C. Kittel, Introduction to Solid State Physics (John Wi- 
ley and Sons, Inc., New York, NY USA, 1996), 7th ed. 

G. Feher, Phys. Rev. 114, 1219 (1959). 

S. Maekawa and N. Kinoshita, J. Phys. Soc. Japan 20, 
1447 (1965). 

Y. Ochiai and E. Matsuura, Phys. Stat. Sol. (a) 38 
(1976). 

G. Wagoner, Phys. Rev. 118, 647 (1960). 

M. Nechtschein et al, Phys. Rev. B 27, 61 (1983). 

K. Matsubara, T. Tsuzuku, and K. Sugihara, Phys. Rev. 

B 44, 11845 (1991). 

K. Mizoguchi et al, Phys. Rev. B 51, 8864 (1995). 

A. J. Epstein et al, Syn. Metals 21, 63 (1987). 

Q. Li, L. Cruz, and P. Phillips, Phys. Rev. B 47, 1840 

(1993). 

H. Ue and S. Maekawa, Phys. Rev. B 3, 4232 (1971). 
M. N. Alexander and D. F. Holcomb, Rev. Mod. Phys. 
40, 815 (1968). 

J. D. Quirt and J. R. Marko, Phys. Rev. B 5, 1716 (1972). 

C. Gomez-Navarro et al, Nat. Mat. 4, 534 (2005). 

F. Demichelis et al, Diam Relat. Mater. 3, 844 (1994). 



